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Abstract 

We analyse numerically manifestations of the radiative amplitude zero (RAZ) ef- 
fect in single leptoquark production associated with hard photon emission. We present 
some quantitative conclusions on the possibility to distinguish leptoquark charges 
produced in ep collisions taking account of three-body final state subprocesses and of 
proton structure functions. Applying this method to HERA and possible LEP+LHC 
experiments we show that the RAZ analysis can serve as a tool to determine the 
leptoquark electric charge up to large leptoquark masses. 



(To be submitted to Physics Letters) 



1 Introduction. 



In paper M we have proposed a method for distinguishing the charges of leptoquarks]] (LQ) produced in 
the reaction 

e ± +p-^~/ + LQ + X. 

The idea is based on the radiative amplitude zero (RAZ) effect, i.e. the absence of photon emission at 
some directions depending on the leptoquark electric charges. In this paper we give a detailed numerical 
analysis of this method applying it to HERA and possible LEP+LHC experiments. 

In Ref. || we have derived the analytical formulas for the unpolarized cross sections of the hard 
subprocesses 

e ± + q 7 + LQ (1) 

We have shown that they are proportional^ to a factor (Q e u — Q q t) 2 . So, at some photon emission angles 
where this factor vanishes the cross sections have exact zeros whose positions are completely determined 
by the electron (positron) and quark charges, Q e and Q q , respectively. 

Of course in reality the leptoquarks decay into two fermions. Considering the electron decay channel 
we have got the following clear signature: the hard photon, the electron (or positron) and a hard (quark) 
jet with complete kinematical reconstruction of the subprocess: 

e ± +q -> 7 + e ± +q (2) 

In this paper we analyse numerically the manifestations of the RAZ effect for the case of this three- 
body final state subprocess. We make some quantitative conclusions on a possibility to distinguish the 
charges of leptoquarks produced in ep collisions taking account of the proton structure functions. 

Our calculations were made with the help of the CompHEP package p| created for automatic calcu- 
lation of cross sections at tree level in the Standard Model and beyond. For the numerical integration 
over three-particle phase space and for event generation we have used the BASES package W. 



2 Calculation framework 

The amplitude of subprocess e~ + d — > e~ + d + 7 is represented by the Feynman diagrams displayed in 
Figs. 1, 2. In Fig. 1 the diagrams corresponding to the signal LQ production process (2) are represented. 
The diagrams of Fig. 2 correspond to the SM contribution: deep inelastic scattering associated with hard 
photon emission. 

Using the CompHEP package we have calculated analytically the unpolarized squared matrix elements 
for subprocesses (§). But in contrast to our previous paper jlj we have introduced finite widths into the 
LQ propagators. Then we have evaluated numerically the contribution of separate (quark) constituents 
to the integrated cross sections by the convolution of the subprocess cross section with the corresponding 
parton distribution function: 

a{s) = J dxq(x,Q 2 ) J d$ 3 — e cu4s (£ 7 ,tf 7 ,...). (3) 

Here d&s is the element of 3-body phase space; a is the cross section of subprocess (^]); s is the squared 
CMS energy of the electron-proton system; the squared CMS energy of the hard subprocess equals s = xs; 
the quark distribution function is denoted by q(x,Q 2 ) and the 4-momentum transfer scale is taken to be 
Q = V§. The photon emission angle is denoted by i? 7 . The direction i9 7 = is along the proton beam. 
As the lower bound x m i n we use the value M 2 /s to take into account the quasi-resonant peak, where M 
is the LQ mass. The function Q cuts (Ej, . . .) corresponds to the kinematical cuts which we impose to 
get a realistic distribution and values of cross sections. In particular this function has to include cuts 
on the photon energy and on the angles between the photon and all charged particles in (2) to remove 
infrared and collinear divergences. It is also necessary to apply cuts for a reliable separation of the hard 
jet. 



^^The problem of identification of LQ is discussed now in the literature, see recent publication [g| where 
the determination of the leptoquark properties in polarized e7 collisions is considered. 
2 For the vector leptoquark this is correct only in the case of Yang-Mills coupling. 
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All the fcrmion masses are put equal to zero in our calculations. 

The numerical analysis of the corresponding cross sections is given for two cases: 

• HERA: y/s — 296 GeV, electron beam energy E e — 26.7 GeV, integrated luminosity L = lOOpb -1 ; 

• LEP+LHC: V« = 1740 GeV, £ e = 100 GeV, annual integrated luminosity L = 1000 pb -1 . 

One of the general restrictions on the LQ-fermion interaction is the chirality of the fermion-LQ coupling 
(see Q and references therein). Therefore we consider only either couplings with left-handed leptons or 
with right-handed ones. In both cases we use in the numerical analysis an "electroweak" value for this 
constant, A = 0.3. 

As to LQ mass we analyse cross sections in the following ranges 

• HERA: 150 GeV < M < 290 GeV; 

• LEP+LHC: 200 GeV < M < 1500 GeV. 

For the parton densities we used the parametrizations STEQ2p |(| and MRS- A 0, which take account 
of recent HERA data. Both parametrizations gave the same results within calculation error. 
Let us describe in more detail the cuts which we apply in the laboratory frame. 

For a proper background analysis we introduce kinematical cuts for this reaction in such a way as 
to suppress contributions from the Standard Model Feynman diagrams (see Fig. 2). Recall that for the 
problem under consideration the LQ mass is supposed to be already known from resonant LQ production. 
So we can introduce a narrow cut on the invariant mass of the outgoing electron (or positron) and the 
(quark) jet: 

M ut + P° q ut ? ~ M 2 \ < 6MT, (4) 

where T is the LQ width. We have calculated T at tree level for each value of the LQ mass. 

We introduce also a cut on the angle between the outgoing electron (positron) and the quark: d eq > 
10°. In the case of LQ decay the products are moving back-to-back in the LQ rest frame, so this cut 
suppresses the standard model contributions and does not change the LQ signal. 

Then wc have to introduce cuts on the energies and emission angles of the outgoing electron (positron) 
and quark to exclude the forward and backward cones and various soft contributions: 

E e , E q > 10 GeV, 10° <& e ,ti g < 170°. (5) 

We also apply cuts on the photon energy and emission angle: E 1 > 1 GeV, 10° < ?9 7 < 170°. 

Finally we introduce cuts on the angles between the outgoing photon and the electron (positron) and 

quark: i? 7e > 10°, d lq > 10°. These cuts remove the collinear divergences. 

For this set of cuts and with x m i n — (200 GeV) 2 /s the cross sections of reaction (Q) calculated in the 

Standard Model are: 

HERA ct sm ~ 0.5 fb; 

LEP+LHC crg M - 40 fb. 

So, for our analysis we have an SM background of less than 1 event for HERA (assuming an integrated 
luminosity of 100pb _1 ) and less than 50 events per year for LEP+LHC. For larger LQ masses from the 
corresponding ranges (see above) we introduce larger cuts x m i n , so the SM background is smaller. 

The final results presented in this paper were calculated for reaction (Q) with a complete set of tree 
level diagrams, i.e. taking into account the SM diagrams, the contributions of the signal reaction ([!]) and 
the SM-LQ interference terms. 



3 Numerical results and RAZ analysis 

In this section we present results of the numerical analysis for the left-handed sector of LQ-fermion 
interactions. For the right-handed sector the results for all channels are the same as for the corresponding 
processes in the left-handed sector due to identical analytical formulae For the leptoquarks we use 
the notation of || (see also the detailed table of LQ quantum numbers in Jl[]). 

First we note that the cross sections of interest are of order 1 pb for HERA and 10 pb for LEP+LHC. 
These values arc obtained for reaction (0) with the kinematical cuts discussed in section and with 
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M = 200 GeV for HERA and M = 300 GcV for LEP+LHC with an electroweak value of the fermion-LQ 
coupling constant A = 0.3. So, there are about one-hundred events at HERA and more than ten thousand 
events at LEP+LHC. This must be compared with the SM background, which is less than 1 event at 
HERA and 50 events per year at LEP+LHC. 

A typical example of a distribution in x is shown in Fig. 3, for the S^ 1 leptoquark with a — 0.82 pb 
for HERA at M = 200 GeV and a = 13 pb for LEP+LHC at M = 300 GeV. One can see that there is a 
sharp peak near x m i n — M 2 /s. 

The cross sections decrease with increasing leptoquark mass. This is shown in Fig. 4 for HERA and 
in Fig. 5 for LEP+LHC. 

In Fig. 6 we present the angular distributions of hard photons for the case of scalar leptoquarks and 
u and d valence quarks. Again we present these distributions at Mlq = 200 GeV for HERA and at 
Mlq = 300 GcV for LEP+LHC. 

For HERA we took 16 bins, 10° per bin. Due to the small statistics (even for an integrated luminosity 
of 100 pb -1 ) a finer binning does not allow us to recognize the RAZ effect. For LEP+LHC we took 40 

bins, 4° per bin. We see explicit RAZ for S^ 1 and R 2 2 leptoquarks at different values of RAZ angles in 
good agreement with Ref. [Q where the RAZ angles were estimated analytically. For the e~u channel 
(leptoquarks Si, S®) this angular distribution has no radiative amplitude zero^ and the distribution in 
the e + d case is similar. 

The angular distributions in for vector leptoquarks with Yang-Mills couplings and valence quarks 
in the initial state have a behaviour similar to the corresponding cases for scalar LQ; they are therefore 
not shown. 

Also for channels with sea quarks the distributions in the photon emission angle are similar but with 
a significantly smaller level of cross sections (see Fig. 5 ) . Note that for HERA there is no statistics for 
any of the processes in the channels with sea quarks. 

Considering these distributions we see that for HERA the RAZ analysis seems to be statistically 
unsupported in the case of leptoquark masses greater than 200 GeV (for A = 0.3). However from the 
HERA data analysis corresponding to an integrated luminosity of ~ 425 nb" 1 |g[ for the i?2 and i?2 
scalar LQ and for the Ui, U\ and U3 vector LQ the established mass limit is smaller than 200 GeV (for 
A = 0.3). In this small mass range between the experimental limit and 200 GeV the determination of the 
LQ electric charge is still possible by the proposed method if this new boson is discovered at HERA. 

For LEP+LHC there is enough statistics for a reliable RAZ analysis up to large leptoquark masses. 
Consider, as a criterion, 10 events per 4° bin, in which case we can say with confidence whether RAZ 
is present or absent. There are 40 bins in the interval 10° < ?9 7 < 170°, so the lower limit for the cross 
section is 0.4 pb to be certain that RAZ is statistically visible. From the leptoquark mass dependence of 
the cross sections (see Fig. 4 and Fig. 5) we derive the upper limits for leptoquark masses with visible 
RAZ. These mass limits are shown in Tabic [TJ. 

We see that for LEP+LHC the RAZ analysis, and indeed a measurement of the leptoquark charge, can 
be made in the e~p and e + p channels separately up to M ~ 500 — 600 GeV for all types of leptoquarks. 
For higher masses, up to 850 — 900 GeV, all types of leptoquarks could also be identified but only some 
of than in the e~p channel and others in the e + p channel. The upper limit for the RAZ analysis is about 
M ~ 1 TeV, see details in Table |l|. For the vector case the situation is slightly better for U3 and U\ 
leptoquarks. Here the RAZ effect can be reliably determined up to M ~ 1200 GeV. 

Note that the pairs (Si, S®) and (U\, U®) (in the left-handed sector) will remain unresolved by this 
method, because of their equal charges and third components of isospin. 

Finally we note that our analysis was made for the reaction with the electron decay mode of LQ. 
So only leptoquarks with the third components of isospin equal to —1 and — ^ from the corresponding 
isotriplets and isodoublets were considered in the analysis (with the only exception of R 2 and V 2 in the 
right-handed sector). To analyse leptoquarks with complementary isospin projections from the same 
multiplct it is necessary to investigate the neutrino decay mode of LQ. However the signature of these 
processes includes missing energy due to the neutrino emission. So one has to carry out the analysis 
differently in this case. This analysis will be presented elsewhere. 

3 The drop in the bin 130° < i? 7 < 140° in this channel is a statistical fluctuation. 
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Conclusions 



Our principal conclusion is that for some types of leptoquarks (in the left-handed sector: S 3 and R 2 2 

scalars, V 2 2 and U 3 vectors; in the right-handed sector: S\ and R 2 2 scalars, V 2 2 and U\ vectors) the 
radiative amplitude zero is present in the distribution in the photon emission angle. This effect can be 

detected at LEP+LHC up to large values of leptoquark masses: 900 GeV for S^ 1 and V 2 2 , 1 TeV for 

#2^ and 1.2 TeV for U^ 1 and U x . 

The RAZ analysis (whether the radiative amplitude zero is present or not) gives us a possibility to 
determine the leptoquark electric charge (and, as a result, its other quantum numbers). This method is 
available up to rather large leptoquark masses (up to 1 TeV at LEP+LHC). Even at HERA there still 
exists a small mass region (near 200 GeV) for some leptoquarks where the RAZ effect can be observed. 

For vector leptoquarks the exact RAZ effect is present only in the case of Yang-Mills type of photon- 
LQ coupling. So we have got a fairly sensitive tool to measure the anomalous magnetic moment n in the 
case of RAZ. A detailed analysis of this possibility will be presented in a future publication. 

All our numerical results were obtained for the electroweak value of the fermion-LQ coupling constant, 
A = 0.3. For other values of A numerical estimates can be obtained by rescaling our results (cross sections 
have A 2 as a factor). 
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Figures captions 

Figure 1: Feynman diagrams with LQ signal. 

Figure 2: Feynman diagrams for subprocesses e~ + q — > e~ + q + 7 in SM. 
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Figure 3: Distributions in x for e +d^j + e +d with the S 3 1 leptoquark contribution. These distributions 
are obtained at M = 200 GeV for HERA (a = 0.82 pb) and at M = 300 GeV for LEP+LHC 
(a = 13 pb) with kinematical cuts applied. 

Figure 4: Cross sections in dependence on the LQ mass at HERA. 

Figure 5: Cross sections in dependence on the LQ mass at LEP+LHC. 

Figure 6: Distributions in the photon emission angle $ 7 for valence quarks (M = 200 GeV for HERA and 
M = 300 GeV for LEP+LHC) in the case of scalar LQ. 



5 



£q channel 


e u 


e~d 


e~u 


e~d 


e + u 


e+d 


e+u 


e + d 


scalar LQ 


L 




Si, S3 


s; 1 


i 

R 2 2 




1 

i? 2 2 


i? 2 2 


Si,S$ 


5T 1 






R 


Si 


Si 


_ i_ 
R<2 2 


1 


_ 1 


1 


Si 


Si 


max M [GeVJ 


900 


900 


540 


610 


980 


830 


670 


500 


vector LQ 


L 






v;\ 


ui 1 


Ui,U$ 


Us 1 






v 2 \ 






R 




VP 


Ui 


Ui 


Ui 








max M 


GeV 




1140 


910 


610 


690 


1180 


830 


660 


660 



Table 1: LEP+LHC. The RAZ upper limits for LQ masses (in the positron channels anti-LQ are pro- 
duced) . L-row corresponds to the interaction of LQ with left-handed leptons and R-row to right-handed 
ones. 
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